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Abstract 

By means of the mixed representation RPA based on the Skyrme-Hartree-Fock 
mean field, we investigate low-frequency octupole excitations built on the superde¬ 
formed (SD) states in the N = Z nuclei around ^^Ca and the neutron-rich Sulfur 
isotopes. The RPA calculation is carried out fully self-consistently on the three- 
dimensional Cartesian mesh in a box, and yields a number of low-frequency octupole 
vibrations built on the SD states in ^^S, ^®Ar, ^^Ca and ^^Ti. In particular, a strongly 
collective = 1“ octupole vibration is suggested to appear on top of the SD state 
in ^°Ca. For neutron drip line, we find that the low-lying state 

created by the excitation of a single neutron from a loosely bound low D level to a 
high D resonance level acquires an extremely strong octupole transition strength due 
to the spatially very extended structure of the particle-hole wave functions. 


1 Introduction 

Nowadays, more than two hundreds superdeformed (SD) bands are identihed in various 
mass (A=60, 80, 130, 150, 190) regions [1-6]. Every SD regions have their own characteris¬ 
tics so that we can significantly enlarge and deepen our understanding of nuclear structure 
by systematically investigating similarities and differences among the SD bands in different 
mass region. The SD shell structure is significantly different from that of normal defor¬ 
mation. Namely, each major shell at the SD shape consists of about equal numbers of 
positive- and negative-parity levels. This is a favourable situation for the appearance of 
negative-parity collective modes. In fact, various mean-held calculations [7-9] and quasi¬ 
particle random phase approximation (RPA) [10,11] on the basis of the rotating mean held 
(cranked shell model) indicated that SD nuclei are very soft against both the axial and non- 
axial octupole deformations. Thus, low-frequency octupole vibrations have been predicted 
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to appear near the SD yrast lines [12], and recently discovered in experiments for heavy SD 
nuclei in the Hg-Pb region [13], and also in ^^^Dy [14], 

In recent years, the SD bands have been discovered also in the ‘^‘’Ca region: In ^®Ar 
the SD band has been identihed up to its termination at = 16^ [15-17]. The SD band 
in the spherical doubly magic nucleus ^'^Ca is built on the 8p-8h excited O’*" state at 5.213 
MeV [18,19]. The SD shell gap at iV = 20 (Z = 20) is associated with the neutron 
(proton) 4p-4h excitation from below the iV = 20 (Z = 20) spherical closed shell to the 
fj /2 shell. The rotational band built on the excited O’*" state at 1.905 MeV in ^"^Ti may 
also be regarded as belonging to a family of the SD band [20]. The fact that rotational 
bands built on excited O’*" states are systematically observed is a unique feature of the SD 
bands in the ‘^‘^Ca region, considering that the low angular momentum portions of the SD 
bands in heavy nuclei are unknown in almost all cases (except the fission isomers). We have 
conhrmed that the symmetry-unrestricted Skyrme-Hartree-Fock (SHF) calculation indeed 
yields the SD local minima corresponding to these experimental data [21]. It should be 
stressed that, in spite of the remarkable progress in experiment, the doubly magic SD band 
in ^^S associated with the SD magic number N=Z=16, which has been anticipated quite a 
long time [22-31], has not yet been observed and remains as a great challenge. 

If we believe that Z=16 is a good SD magic number, the existence of the SD shell gap 
at TV = 20, revealed by the discovery of the SD band in ^°Ca [18], suggests that another 
SD band would appear in the neutron-rich nucleus ^®S. Furthermore, combining with the 
fact that the SD bands have been observed in [32,33], we can expect that Sulfur 

isotopes around ^®S, which are situated close to the neutron-drip line [34,35], constitute 
a new SD region associated with the SD shell gaps at Z = 16 and iV ~ 32. In fact, the 
symmetry-unrestricted SHF calculation [36] yields the SD local minima for as well as 

32S. 

The investigation of low-frequency octupole vibrations built on the SD states in the 
yl=30-50 region possesses some new features that are absent in the study of heavy SD 
nuclei. For the N = Z nuclei in the ^°Ca region, it may be possible to observe in experiment 
such collective modes built on the known SD O’*" states. Moreover, because the proton and 
neutron shell structures are essentially the same, we can expect that strong coherence takes 
place between the proton and neutron excitations and brings about an enhanced collectivity 
of these modes. Concerning the anticipated new SD region around ^®S, we will encounter 
an essentially new situation: Because these nuclei are situated close to the neutron drip 
line, there is almost no bound state above the Fermi surface and nucleons may excite into 
the continuum states even in the lowest excited states. 

The study of soft collective modes unique to unstable nuclei close to the neutron drip 
line is one of the current major subjects in nuclear structure physics. For this purpose, 
the continuum RPA [37, 38] based on the Green’s function method correctly treating the 
continuum boundary condition has been widely used (see, e.g., [39,40]). Recently, this 
method has been extended [41-43] to include the pairing correlations on the basis of the 
Hartree-Fock-Bogoluibov theory. It is called the continuum quasiparticle-RPA (QRPA). 
Furthermore, a new technique of solving the QRPA equations using the canonical basis has 
been developed [44]. The canonical basis has also been used in the QRPA calculation based 
on the relativistic mean held scheme [45,46]. A serious limitation of these works is that, 
because the construction of the Green’s function for deformed mean helds is difficult, they 
are restricted to spherical nuclei. Quite recently, however, this limitation has been overcome 
by Nakatsukasa and Yabana [47]; they have proposed an iterative method of constructing 
the response functions for deformed systems with the proper boundary condition in the 
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three-dimensional (3D) coordinate space. They have also proposed a feasible method of 
treating the continuum boundary condition by means of the real-space time-dependent HF 
approach with the absorbing boundary condition. 

In this paper, we investigate low-frequency octupole excitations built on the SD states in 
neutron rich Sulfur isotopes as well as those in the ^°Ca region. For this purpose, we employ 
the mixed representation RPA [48-50] based on the SHF mean held. In this RPA scheme, 
the particles above the Fermi surface are described in the coordinate representation, while 
the holes are represented in the HF single-particle basis. The RPA calculation is carried 
out fully self-consistently on the 3D Cartesian mesh in a box. The major advantage of 
this approach is that it is a fully self-consistent scheme in the sense that the same effective 
interaction is used in both the mean held and the RPA calculations; i.e., all terms of the 
Skyrme interaction contributing to the RPA equations of motion are taken into account. 
It should be stressed here that such a fully self-consistent calculation using the Skyrme- 
type interaction is difficult in the Green’s function approach, so that, usually, the residual 
particle-hole interactions associated with the spin-orbit and Coulomb interactions, as well as 
those associated with the time-odd components (related to the spin and current densities) of 
the mean held, are ignored in the continuum RPA based on the SHF mean held. Another 
important merit of our approach is that, thanks to the use of the 3D Cartesian mesh 
representation, we can treat strongly deformed nuclei on the same footing as spherical 
nuclei. Furthermore, because the particles are described in the coordinate representation, 
we do not need to introduce an upper cut oh for their energies. Quite recently, Imagawa 
and Hashimoto constructed a new computer code that carries out the fully self-consistent 
RPA calculation in the mixed representation on the basis of the SHF mean held, and they 
carefully tested its numerical accuracy [51-53]. The numerical calculations of this paper 
were carried out using a rehned version this code. 

This paper is organized as follows: In Section 2, a brief account of the self-consistent 
RPA calculation using the mixed representation is given. In Section 3, we present and 
discuss the results of numerical calculation for low-frequency octupole vibrations built on 
the SD states m. N = Z nuclei, ^^S, ^®Ar, ^°Ca, and "^^Ti. In Section 4, the results for 
neutron-rich Sulfur isotopes, ^®S, ^®S, and ^°S, are presented and discussed. The main 
results of this paper are summarized in Section 5. 

A preliminary version of this work was reported in Refs. [54,55]. 


2 RPA calculation using the mixed representation 

2.1 Basic formulae 

The RPA equations in the mixed representation are easily derived either by means of the 
lineralized equations of motion approach [48] or the small-amplitude approximation of the 
time-dependent HF theory [49]. Here, we recapitulate the basic formula following the former 
approach. 

The creation and annihilation operators of nucleon, '^^(a;) and 'ip{x), are divided into 
those for the bound states and those for the continuum states: 

COO 

i’Kx) =Y.‘h’a{x)ci+ depl{x)c\, ( 1 ) 

Q io 

coo 

=^iPa{x)Cc+ j <leipe{x)Ce, ( 2 ) 

V -fo 
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where x denotes a set of space, spin and isospin coordinates, i.e., x = {r, a, r}. The symbol 
a denotes the HF bound states with wave functions (Pa(x), while the energy e specihes 
the continuum states with wave functions (ps(x). With the use of the particle and hole 
concepts, we then divide the nucleon operators according to the occupation number da in 
the HF ground state {9a = 0 for unoccupied states, 9a = I for occupied states). We consider 
even-even nuclei and assume that the single-particle states a and their time reversed states 
are doubly degenerated. Obviously, the continuum states belong to the particle space. On 
the other hand, for the bound states, the operators (c)^, c^) are divided into the particle 
operators {a}^, am) and the hole operators (6|, bi), 

4 = (1 - (3) 

c« = (1 - 9a)am + 9abl, (4) 

where the indices m and i are used in place of a to distinguish the particle and hole states. 
Thus, the particle creation and annihilation operators, a^{x) and a{x), at the coordinate x 
are written in terms of the above {a}^, am) and the integration with respect to the positive 
energy e specifying the continuum states; 

m 

= ^‘^m{x)am + / de(pe{x)a^ 

m 

In the mixed representation, the coordinate representation is used for particles, while the 
HF basis specihed by the discrete index i are used for holes. Thus, the following expressions 
of the nucleon operators are convenient: 


'^\x) = a\x) + Y.ip*{x)bi, 

% 

(7) 

^{x) = a{x) + '^^i{x)b\. 

(8) 


i 


Note that 

{a{x),a\x')} = P{x,x') = 5{x,x) - Q{x,x') = 5{x,x) - Y^(pi{x)(p*{x), (9) 

i 

where P{x, x') and Q{x, x') are the projectors onto the particle and hole spaces, respectively. 

Using an abbreviation, / d?", the RPA phonon creation operators in the mixed 

representation are written as 

- Y^^{x)bia{x)} . ( 10 ) 

i X 

Then, the RPA eigenvalue equations determining the eigen-energies, hu\, and the forward 
and backward amplitudes, X^{x) and Y^^{x) are given in the following form: 

= huJxXi{x) ( 11 ) 

j x' 

J2T.lB’„{x,x')Xpx') + A-^(x,x')Ypx')] = -huixYpx), (12) 

j x' 


(5) 

( 6 ) 
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where 


Ay(x,x') 

= P{x,x")Aij{x",x"')P{x"\x'), 

x" ^x’’’ 

( 13 ) 

Bij{x,x') 

= Y P{x,x'')B^{x!',x!'')Q{x!'',x), 

x''x"' 

( 14 ) 


and 


Aij{x,x') 


Bij{x,x') 


[hnpix, x') - eiS(x, x')] 5ij + ^ (pi{xA^*j{x" 


d^E[p] 


dp{x”, x)dp{x', x'") ’ 




d‘^E[p] 


dp{x", x)dp{x'", x') 


( 15 ) 

( 16 ) 


The above RPA equations can be recast into the following form for the linear combinations, 
4p\x) = Xpx)±Yp(xy. 


^ \hYa(x,x') - eiS(x,x')]<l>P^(x') + P(x,x')h'-^^px',x'')(f,(x") = hujxil>[ ^(x) , (17) 

x' x'x" 

^ [h]iY{x,x') - ei5{x,x')] ^ P{x,x')h‘^^^l{x',x")^i{x") = hujx(t)t^^{x) . (18) 




Here, the RPA Hamiltonian is given by 


K^a[PX)- p p [X^X 


(19) 


where 


p^^P(x, x') = Y, ± ‘p‘'P’(x')f>i{x) 


( 20 ) 


represents the transition density. The above equations are obtained also by making the 
small-amplitude approximation for the equations of motion of the time-dependent HP mean 
held [49]. The explicit expression of hj(^^(a;, x') is given in Appendix. The transition matrix 
elements for the octupole operators Qsic are given by 

(OlQsxIA) = Y.Q^K{r)Y^ [Xt{x)p*{x) + pi{x)Y^^{x) . 

X i 


2.2 Details of numerical calculation 

We solve the SHF equation on the 3D Cartesian mesh assuming the rehection symmetry 
about the (x, y)-, (y, z)-, and {z, a;)-planes. The derivative operators are treated by means 
of the Lagrange mesh method [56], except for the derivative of the Coulomb potential 
for which the hnite difference method with the 9 points formula is used. Based on the 
numerical test explained below, we adopt the mesh spacing Ax = 0.6 fm. Because we treat 
the superdeformed nuclei, we use the rectangular box with 15 and 25 mesh points in the 
positive directions of the minor and major axes, respectively, i.e., we take the mesh points 
at X = 0.3, 0.9,..., 8.7 fm, y = 0.3, 0.9,..., 8.7 fm and x = 0.3, 0.9,..., 14.7 fm (the major axis 
is called the x-axis). Considering the reflection symmetries, the total size of the box is thus 
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17.4 X 17.4 X 29.4 fm^. For the effective interaction, the standard versions of the Skyrme 
interaction, Sill [57], SkM* [58] and SLy4 [59], are used. 

We then solve the RPA equations, (17) and (18), using the SHF solutions obtained 
above. With the mesh discretization, Eqs. (17) and (18) become the matrix eigenvalue 
problem with the 2 x W x Np dimension, where Nh and Np denote the numbers of the hole 
states and the mesh points, respectively. Because we assume the reflection symmetry about 
the three planes, the RPA eigen-value equations are separated into four blocks specihed by 
the parity quantum numbers. It is easily ascertained that the negative parity operators 
are classihed into the four groups, (—,-|-,-|-), (-|-,—,-|-), (-|-,-|-,—), and (—,—,—), where 
the hrst, the second, and the third signs indicate the parity with respect to the reflection 
about the {y, z)-, {z, x)-, and (x, |/)-planes, respectively. For instance, the octupole operator 
Qso = belongs to the (-I-, -7, —) sector. On the other hand, the octupole operators 

with iP ^ 0, gg = r3 {Y,k ± Y,,-k) /v^, are classified according to the parity quantum 
numbers as shown in Table 1. We calculate the eigenvalues and eigenvectors of the RPA 
matrix by means of the conjugate gradient method in a version developed in Refs. [51-53]. 

It should be noted here that the translational symmetry is broken in the coordinate 
mesh representation so that the eigenvalues of the RPA corresponding to the spurious 
center of mass motion do not necessarily become zero. It is possible that the eigenvalue 
huj\ corresponding to the spurious mode become pure imaginary. Moreover, in deformed 
nuclei, except for the (—, —, —) sector, these spurious modes can mix with the octupole 
modes having the same symmetry properties (see Table 1 for the classification of the dipole 
operators corresponding to the center of mass motion). It is therefore very important to 
carefully check the accuracy of numerical calculation whether or not the spurious center 
of motions are decoupled from the physical excitations in a good approximation. In the 
numerical calculation, we manipulate this problem with the following procedure. Combining 
Eqs. (11) and (12), the RPA eigenvalue problem is transformed into the following form: 

jk x' ^x" 

YBij{x,x)Ajk{x^- Bij{x,x)Bji,{x\x^')\^^^^^{x') = ( 21 ) 

This form is convenient for applying the conjugate gradient method which is valid only for 
real eigenvalues, because, when the eigenvalue hc^A of the spurious mode take an imaginary 
value, we can easily identify it as a negative {huj\Y solution of this equation. The result 
of such an accuracy test of numerical calculation is presented in Fig. 1. In this hgure, 
low-energy solutions of the above RPA equation for negative-parity excitations on the SD 
ground state in ^°Ca are plotted as functions of the mesh spacing Ax used in the SHF- 
RPA calculation with the Sill interaction. We see that a good convergence is attained 
at about Ax = 0.6 fm, although excitation energies of the spurious center of mass modes 
fluctuate in the region of Ax >0.7 fm. It should be noted that the excitation energies of 
the physical excitation modes and their transition strengths almost converge for Ax <0.7 
fm. Thus, we adopt Ax = 0.6 fm in the numerical calculations presented below. We note 
that the use of this mesh spacing corresponds the introduction of an effective energy cut-off 
.Ecut-oflf = /2m){ tt/A x)'^ ~ 500 MeV. This value is sufficiently large. Accordingly, it may 
be allowed to say that our calculation is, in practice, cut-off free. Certainly, it is desirable to 
use a smaller Ax and a larger box for better numerical accuracy, but it is difficult to do so 
with the available computing power at the present time. This limitation will be overcome 
in the near future. 
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In presenting the results of the SHF-RPA calculation in the subsequent sections, we label 
the SHF single-particle levels in terms of the asymptotic quantum numbers [iVosc, A, H] 
for the modihed oscillator potential. It should be noted that this labeling is used solely for 
convenience of presentation, and do not always mean that they are good quantum numbers. 
The labeling of the SHF single-particle states is done by examining their properties, like 
reflection symmetries and expectation values of the angular momentum, and by comparison 
with the single-particle level scheme obtained in Ref. [60], where the deformed Woods-Saxon 
potential simulating the SHF potential was constructed and properties of both the bound 
states and the discretized continuum states for this potential are analysed in detail. In this 
reference, the RPA calculations were also carried out for the octupole excitations on the 
SD states under consideration, using the conventional matrix formulation with these single¬ 
particle levels in a truncated model space and a density dependent contact interaction. 
The results of this simple RPA calculation qualitatively agree with the results of the mixed 
representation RPA calculation presented below. Thus, we have used the correspondence 
between the two calculations for checking the numerical calculation. As we shall discuss 
later, this comparison was especially useful for the purpose of distinguishing the resonant 
levels from the non-resonant discretized continuum levels. 

3 Results of calculation for the ^“^Ca region 

Figure 2 displays the transition strength distributions obtained by the SHF-RPA calculation 
with the SHI interaction for low-frequency octupole excitations built on the SD states in 
^^S, ^®Ar, and ^°Ca. The RPA transition stengths are here dehned as the squared 

matrix elements of the isoscalar octupole operators between the SD excited state |A) and 
the SD ground state |0) , = \ {^\Q'Ik^Q'^k\'^)\^■ hgure, the RPA transition 

strengths for the electric octupole operators, B{E3) = |(0|(53if|A)p, are also indicated for 
levels possessing strong collectivity. Here, and denote the proton and neutron 
octupole operators. It should be noted that these quantities represent the intrinsic matrix 
elements squared. Accordingly, in order to evaluate the reduced transition probabilities 
in the laboratory frame, appropriate Clebsch-Gordan (CG) coefficients squared and other 
kinematical factors should be multiplied to them [61]. Specifically, for the transition from 
the 3“ member of the K = 0 {K ^ 0) band to the SD O’*" band head, the GG coefficient is 
unity and the kinematical factor is 1/7 (2/7). The SHF calculation for ^°Ga in fact yields 
a SD mean held with small triaxially of 7=7°, while the SD solutions for ^^S and ^®Ar 
are axially symmetric. Nevertheless, for simplicity, we here present the result of the RPA 
calculation on the basis of the SHF mean held obtained under the constraint 7=0°. The 
ehect of the small triaxial deformation on these excitation modes will be discussed in the 
end of this section. 

For ^^S, we obtain a low-frequency collective K'^ = 2“ mode at 2.6 MeV above the 
SD band head. It possesses a large isoscalar octupole transition strength of 30 W.u. (1 
W.u. ~ 61 fm® for ^^S). The unperturbed strengths in the 3.5 MeV region, seen in Fig. 2, 
are associated with the proton and neutron excitations from the [211|] level to the [3211] 
level. The RPA transition strength is greatly enhanced in comparison with the unperturbed 
strengths and the RPA eigen-energy is signihcantly shifted down from the unperturbed 
particle-hole excitation energies, so that this K'^ = T~ mode can be regarded as a collective 
octupole vibration. 

For ®®Ar, we obtain two strongly collective modes; a K'^ = 2“ mode at 3.9 MeV and 
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a K'^ = 1“ mode at 4.0 MeV. They possess large isoscalar octupole transition strengths 
of 27 and 25 W.u., respectively (1 W.u. ~ 77 fm® for ^®Ar). This 2“ mode is similar to 
that in discussed above, except for the following difference: Because the quadrupole 
deformation of the SD state in ^®Ar is smaller {P 2 = 0.51) than that in {P 2 = 0.72), 
the unperturbed particle-hole excitation energies (4.9-5.0 MeV) of protons and neutrons 
from the up-sloping [21l|] level to the down-sloping [321|] level are larger in ^®Ar than in 
Accordingly, the RPA excitation energy of the 2“ mode in ^®Ar is higher than that 
in Concerning the 1“ mode, its excitation energy is considerably shifted down from 
the particle-hole excitation energies (4.6-4.9 MeV) of protons and neutrons from the [202|] 
level to the [3211] level. Considering this together with the enhanced transition strength, 
this 1“ mode can also be regarded as a mode possessing strongly collectivity. In addition 
to the 2“ and 1“ modes mentioned above, we also obtain a 0“ mode and another 1“ mode 
in the 2.8-2.9 MeV region, which are moderately collective. 

For ^‘^Ca, we obtain a low-frequency collective K'^ = 1“ mode at about 1.2 MeV. It 
possesses a large isoscalar octupole transition strengths of 29 W.u. (1 W.u. ~ 95 fm® for 
‘^'’Ca). The unperturbed particle-hole strengths in the 2.2 MeV region are associated with 
the proton and neutron excitations from the [32l|] level to the [200^] level. The fact that 
the RPA transition strength is signihcantly enhanced in comparison with the unperturbed 
strength and the RPA eigen-energy is shifted down from the unperturbed particle-hole 
excitation energies indicates that this K'" = 1“ mode possesses a strong collectivity. We 
also obtained a strongly collective K'^ = 0“ mode possessing quite large strength (82 W.u.) 
at about 3.3 MeV. For this K'^ = 0“ transitions, although small unperturbed strengths 
(about 0.1 W.u.) associated with the particle-hole excitations of protons and neutrons 
from the [32l|] level to the [202|] level are seen at about 4.6 MeV, their strengths are small 
mainly because they strongly violate the asymptotic selection rule for low-energy octupole 
transitions in the SD harmonic-oscillator potential [11]. Thus, this collective K'^ = 0“ 
mode is generated mainly from many particle-hole conhgurations lying in the energy region 
higher than that shown in this hgure. We mention that these collective modes in ^®Ca were 
previously reported by one of authors (H.I.) in Ref. [51]. 

For all the RPA eigen-modes presented in Fig. 2, we can say the followings: 

1) The approximate relation, ^ 4R(i73), holds, indicating that they are collective 

modes generated by coherent superposition of proton and neutron excitations. 

2) The RPA strengths are greatly enhanced from the unperturbed strengths of individual 
particle-hole excitations. This indicates that these RPA eigen-modes are generated by 
collective contributions of quite a number of particle-hole excitations. 

The role of coherence between the neutron and proton excitations in building up these 
collective enhancements can also be seen by comparing these RPA transition strengths for 
the N = Z nuclei, and ^®Ca, with those for the N ^ Z nucleus presented in the 
succeeding section. 

We also carried out the same kind of calculation using the SkM* [58] and the SLy4 [59] 
interactions. The results for strongly collective excitations of these SHF-RPA calculations 
with different versions of the Skyrme interaction are compared in Fig. 3. Note that the 
calculations for ^®Ca were done with the constraint 7=0° for the mean helds, although 
the SHF calculations yield small triaxial deformations; 4° and 8° with the SkM* and SLy4 
interactions, respectively. For all excitation modes discussed above, we obtained the results 
qualitatively similar to those with the Sill interaction. Namely, our prediction about the 
existence of these collective modes is rather robust, depending little on the choice of the 
Skyrme interaction. 



Next, let us examine the effects of the triaxial deformations of the SD mean field on 
the RPA eigen-modes discussed above. Figure 4 shows the strongly collective excitations in 
^°Ca and "^^Ti, obtained by the RPA calculation taking into account the triaxiality of the SD 
mean fields. Because K is not a good quantum number due to the K mixing, the K values 
denoted on the right-hand sides of individual levels indicate the maximum components, 
and the transition strengths for these K values are indicated beside the arrows in this 
hgure. A clear consequence of the triaxiality is that the K modes split into doublets. 
This qualitative feature is common to the RPA calculations with the Sill, SkM* and SLy4 
interactions. Concerning the ^°Ca case, however, the small triaxilality (7 < 10°) in the HF 
mean held might be easily diminished by, e.g., the zero-point huctuation in the 7 direction 
and/or pairing correlations. Therefore, it is not clear whether this splitting of the 1“ mode 
is important or not. On the other hand, the triaxial deformations for ^"^Ti, obtained in the 
SHF calculations with the Sill and SLy4 interactions, are signihcantly larger than those 
for ^‘^Ca (see the 7 values indicated at the bottom of Fig. 4). In such cases, it may be 
interesting to examine the theoretical calculation against possible experimental signatures 
of the triaxility. For instance, the calculation with the Sill interaction yields a = 0“” 
collective mode at 2.4 MeV and a = 2“” doublet at 3.0 and 3.5 MeV, which possess 
large isoscalar octupole transition strengths of 27, 24 and 16 W.u., respectively (1 W.u. ~ 
115 fm® for ^^Ti). In addition, we also obtain another 0 mode at 5.3 MeV, two 1 doublets 
at 2.5 and 4.2 MeV, and another 2“ doublet at 3.2 MeV, which are moderately collective 
(not shown in Fig. 4 but displayed in Fig. 5). 

The RPA transition strengths are compared with the unperturbed particle-hole strengths 
in Fig. 5. For the A' = 0 transitions, small unperturbed strengths (about 0.15 W.u.) as¬ 
sociated with the particle-hole excitations of protons and neutrons from the [321|] level to 
the [202|] level are seen at about 3.2 MeV. For the K = 2 transitions, the peaks in the 
4.1-4.3 MeV region in the unperturbed strengths are associated with the particle-hole exci¬ 
tations of protons and neutrons from the [200|] level to the [312|] level. These “iP"" = 0“” 
and “iP’^ = 2“” RPA modes are generated by collective superpositions of not only these 
relatively low-lying configurations but also many particle-hole conhgurations lying in the 
higher energy region. Thus, the present RPA calculation suggests that the appearance of 
the 2“ doublet may serve as a good indicator of the triaxial nature of the SD state of ^"^Ti. 
Experimental search for such a doublet on top of the SD band head seems very interesting. 

4 Results of calculation for the neutron-rich Sulfur 
region 

Figure 6 displays the transition strength distributions obtained by the SHF-RPA calculation 
with the Sill interaction for low-frequency octupole excitations built on the SD states in 
36,48,50Let us first discuss the result of calculation for the neutron-rich ^®S. For this nucleus, 
we obtain a iP^ = 1“ mode at 2.6 MeV and a K'^ = 2~ mode at 4.2 MeV. As they possess 
the isoscalar octupole transition strengths of 13 W.u. and 16 W.u., respectively, they are 
moderately collective (1 W.u. ~ 77 fm® for ^®S). It is interesting to compare this result 
with those for and ^®Ca discussed in the previous section. The 1“ mode in this nucleus 
is similar to the 1 “ mode in ^®Ca, but its transition strength is less than half of the latter. 
This reduction of the collectivity is understood in terms of the weakening of the coherence 
between the proton and neutron excitations when one goes from the N = Z nucleus to 
the N ^ Z nucleus. Looking into the microscopic details, we hnd that, for example, the 
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particle-hole excitations of protons and neutrons from the [32l|] level to the [200|] level 
act coherently to produce the K'^ = 1~ collective mode in ^‘^Ca, but this proton excitation 
is absent in the case of because the [32l|] level is unoccupied. Analogous argument 
also applies to the K'" = 2~ mode. This mode is similar to the 2“ mode in but, its 
transition strength is about half of the latter. This is understood again in terms of the 
weakening of the coherence between the proton and neutron excitations. For instance, the 
particle-hole excitations of protons and neutrons from the [21l|] level to the [32l|] level 
act coherently to produce the K'^ = 2“ collective mode in but this neutron excitation 
is absent in the case of because the [3211] level is already occupied. 

Next, let us discuss the result of calculation for and close to the neutron-drip 
line. For we obtain a K'^ = 1“ mode at very low excitation energy, 0.51 MeV. This 
mode possesses a large octupole transition strength of 31 W.u. (1 W.u. ~ 137 fm® for 
^®S). The RPA eigen-energy is shifted down only slightly (0.14 MeV) from the unperturbed 
energy (0.65 MeV) of the neutron excitation from the [43l|] level to the [310|] level. The 
unperturbed transition strength for this single-particle excitation is extremely large, i.e., 
i?((5^^3) = 7.7 W.u. The major cause of this remarkable strength is understood as follows: 
Both the [431|] and [310|] levels are loosely bound with binding energy about 2 MeV. 
Because, as is well known, wave functions of the loosely bound states are signihcantly 
extended outside of the half-density radius, the matrix element of the octupole operator 
between these conhgurations acquires a large contribution from the region outside of the 
nuclear surface. Thus, such a transition strength can become very large in nuclei close to the 
neutron drip. Although the increase of the RPA transition strength from the unperturbed 
strength indicates the presence of collective effects, it may be appropriate to consider that 
the major character of this mode is of single-particle-type. In contrast to the low energy 
peak at 0.51 MeV discussed above, the peaks above the neutron-emission threshold (indi¬ 
cated by arrows in Fig. 6) are considered to be associated with the excitations to discretized 
non-resonant continuum states. We examined that these peaks in the continuum indeed 
moves when we increase the box size in the numerical calculation. Therefore, we cannot 
give dehnite physical meanings to these peaks. 

For we obtain a K'^ = 2“ mode at 2.5 MeV. Like the K'^ = 1“ mode in 
it is essentially a single-particle-type excitation, although it has extremely large octupole 
transition strength, 55 W.u. (1 W.u. ~ 149 fm® for ®®S). Namely, it is created mainly by the 
neutron excitation from the loosely bound [310|] level to the resonance [422|] level. The 
RPA eigen-energy of this mode is shifted down only slightly (0.5 MeV) from the unperturbed 
excitation energy (3.0 MeV) of this particle-hole conhguration. Because the resonance 
[422|] level is situated near the centrifugal barrier top, its wave function is signihcantly 
extended far out side of the half-density radius of this nucleus. Thus, the unperturbed 
octupole transition strength for this particle-hole excitation takes an extremely large value, 
R((5^®3) = 22 W.u. This resonance interpretation of the [422|] level is largely based on 
the analysis made in Ref. [60] where the deformed Woods-Saxon potential simulating our 
SHF potential was constructed and properties of the discretized continuum states for this 
potential are analysed in detail. It was then found that the energy and the root-mean- 
square radius of this level were stable against the variation of the box size. In the same 
paper, this conclusion was conhrmed also by carrying out the eigenphase analysis for this 
level. The appearance of this resonance level is easily understood as due to the rather 
high centrifugal barrier for the relatively large value of the symmetry axis component of 
the angular momentum, H = 5/2. We conhrmed that this 2“ peak at 2.5 MeV does not 
move but the other peaks moves when we increase the box size in the numerical calculation. 
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Thus, like in except for the 2“ peak, the peaks above the the neutron-emission threshold 
are considered to be associated with the excitations to discretized non-resonant continuum 
states. 

We also carried out the same kind of calculation using the SkM* [58] and the SLy4 [59] 
interactions. The results of calculations with different versions of the Skyrme interaction 
are compared in Fig. 7. For the SLy4 interaction, the SD local minimum does not appear 
in so that the RPA calculation was not done for this nucleus. Except for this, these 
Skyrme interactions yield qualitatively similar results for the K'^ = 1“ and K'^ = 2~ modes. 

To accurately describe the spatially extended wave functions, like the loosely bound 
levels and resonance levels in and discussed above, it is certainly desirable to use a 
bigger box in the numerical calculation. Although it is difficult to signihcantly enlarge the 
box size under the present situation of computing power, we plan to manipulate this problem 
by incorporating the adaptive coordinate method into our scheme of numerical calculation. 
We believe, however, that the qualitative features of the present RPA calculation will remain 
valid. We have not obtained well-developed collective octupole modes for and In 
such unstable nuclei close to the neutron-drip line, there is almost no bound state above 
the Fermi surface and neutrons are excited to the continuum states. It seems difficult 
to develop collective correlations among different particle-hole conhgurations under such 
situations. Quite recently, however, one of the authors (M. Y.) pointed out that the pairing 
correlations play an important role for the emergence of low-frequency collective modes in 
such drip-line nuclei [62]. For studying such pairing effects for collective excitations, it will 
be interesting to extend the mixed representation RPA based on the SHF mean held to 
that based on the SHF-Bogoliubov mean held. 


5 Conclusions 

By means of the mixed representation RPA based on the SHF mean held, we have inves¬ 
tigated low-frequency octupole excitations on the SD states in the N = Z nuclei around 
^°Ca and the neutron-rich Sulfur isotopes. The RPA calculation has been carried out on 
the 3D Cartesian mesh in a box, and yielded a number of low-frequency octupole vibra¬ 
tions built on the SD states in ^^S, ^®Ar, ^°Ca and ^^Ti. In particular, a strongly collective 
K'" = 1“ octupole vibration has been suggested to appear on top of the SD state in "‘^Ca. 
For close to the neutron drip line, we have found that the low-lying state created by 
the excitation of a single neutron from a loosely bound low D level to a high D resonance 
level acquires an extremely strong octupole transition strength due to the spatially very 
extended structure of the particle-hole wave functions. 
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A Explicit expressions of the RPA Hamiltonian 

Here, using the notations of Ref. [63], we recapitulate the quantities appearing in the RPA 
Hamiltonian (19). More detailed expressions are given in Ref. [51]. 

The RPA Hamiltonian is written 


Arpa fa'r') = F^'’^ + Y, (»•) A- '5™' 


+ G£'tr)i„. + EG'iU(>’)T;>, (22) 

t3 

AU = -V • + C/g>" - ^ [v ■ + ig’' • v] + C/gg,,,'*.-. (23) 

Gig" = - E W ■ [Cgy V,] + Eg>" -'- Y [v„Bgh + BghvJ e„, (24) 

II ^ III/ 

The quantities, ^ etc. appearing in the above expressions are obtained by the 

second derivative of the time-even and time-odd energy functionals, and with 
respect to the local densities: 
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(±)A dt^ (±)A 
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The local transition densities PttJ^i fhe local transition spin densities the local tran¬ 
sition kinetic energy densities the local transition kinetic spin densities the 

local transition current densities the local transition spin-orbit current tensors dgg 

are given by 


Ptt3 M = Piu {r.r), 
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47 Je 

= 4S\r,r), 


(34) 

474’’) 

<1 

<1 

5 ." 

74 

J?' 

/j^ _ rj^t ^ 

(35) 

TlfP(r] 

= V . V'474’’. r') 

y* _ y*! ^ 

(36) 

774’’) 

= -\(y-v)p\fP(vy) 
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trl 

= -4(v-v')„474’’.’’' 

) ■ (38) 

The local transition spin-orbit 

currents Jf7^ defined by 

j{±)\ _ y- .(±)A 

(39) 


flUUJ 


where 6^,^^ is the Levi-Civita symbol and is the unit vector in the cu-direction. The scalar 
transition densities and the vector transition densities appearing in 

the above equations are dehned by decomposing the transition density matrix x') = 

p(^)^(r(Tr, rV'r') into the spin-isospin channels: 


rW') = 


'r')5„a'5rr' + ’’O ' <^cJa'?>TT’ 

H=-\ tz=-\ 


,7-*3 
(TCT ' rr' 


. (40) 


In the above expressions, the charge-exchange = ±1 components are included for com¬ 
pleteness, although they do not contribute to the excitation modes considered in this paper. 
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Table 1: Classification of octupole and dipole operators according to the parity quantum 
numbers with respect to reflections about the (y, z)-, (z, x)- and (x, |/)-planes. 
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octupole 

dipole 
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Figure 1: Low-energy solutions of the RPA for negative-parity excitations on the SD ground 
state in ^^Ca, plotted as functions of the mesh spacing used in the SHF-|-RPA calculation with 
the Sill interaction. Lower part: Energies of the spurious center of mass modes belonging to the 
(-1-,—and ) sectors are plotted with open diamonds, filled diamonds, and 

open squares, respectively. When their excitation energies, huj\, take imaginary values, the —IcuaI 
values are plotted for convenience of presentation. The filled circles at about 1.0 and 1.7 MeV 
indicate the lowest physical excitations in the (—, -|-, -f) and (-I-, —, -|-) sectors, respectively. Upper 
part: Transition strengths for the lowest physical excitation at about 1.7 MeV in the 
sector, plotted as a function of mesh spacing. 
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Excitation Energy [MeV] neutrons 



Figure 2; Isoscalar octupole transition strengths for K'^ = 0“, 1“ and 2“ excitations on the SD 
states in (a) (b) ^®Ar, and (c) ^°Ca, obtained by the SHF-RPA calculation with the Sill 

interaction. The unit is fm®. For comparison with the RPA strengths, unperturbed particle-hole 
strengths are also shown in the lower panels for each nucleus. There are no significant strengths 
for K'" = 3“ excitations values in this low-energy region, so that they are not shown. Some 
particle-hole excitations of neutrons near the Fermi surface are drawn by arrows with their K'" 
values in the rightmost part for each nucleus. We have similar excitations also for protons. The 
solid, dashed, and dotted lines indicate positive-parity levels, negative-parity levels, and the Fermi 
surface, respectively. The asymptotic quantum numbers [A^n^An] are indicated for pertinent 
levels. 
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Figure 3: Comparison of the SHF-RPA calculations with the Sill, SkM* and SLy4 interactions 
for collective octupole excitations on the SD states in ^^S, ^®Ar, and ^*^Ca. Only collective exci¬ 
tations having greater than 10 W.u. are displayed here (see Fig. 2 for other excitations 

having smaller strengths). The numbers on the right-hand sides of individual levels indicate the 
K quantum numbers. The numbers beside the arrows indicate the intrinsic transition matrix ele¬ 
ments squared of the isoscalar octupole operator, in units of fm®. Those for the electric 

octupole operator, B{E3), are also given in brackets in units of e^fm®. 



Figure 4: Low-frequency collective octupole excitations on the SD states in ‘^'^Ca and ^^Ti, 
obtained by the SHF-RPA calculation taking into account the triaxial deformation of the mean 
field. Only collective excitations having B{Q^^3) greater than 10 W.u. are displayed here. The 
numbers on the right-hand sides of individual levels indicate the approximate K quantum numbers. 
The numbers beside the arrows indicate the intrinsic transition matrix elements squared of the 
isoscalar octupole operator, B{Q^^3), in units of fm®. Those for the electric octupole operator, 
B{E3), are also given in brackets in units of e^fm®. The Skyrme interaction used and the triaxial 
deformation 7 of the mean field are indicated at the bottom. 


19 




Excitation Energy [MeV] 

Figure 5: The same as Fig. 2, but for octupole excitations on the triaxial SD state in ^^Ti. Note 
that, due to the triaxial deformation, K is not a good number for single-particle levels. The 
asymptotic quantum numbers [Nu-zAn] beside the pertinent levels in the rightmost figure merely 
indicate their main components. 
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Figure 6: Isoscalar octupole transition strengths for K'^ = 0“,1“ and 2~ excitations on the 
SD states in (a) (b) and (c) obtained by the SHF-RPA calculation with the Sill 

interaction. The unit is fm®. For comparison with the RPA strengths, unperturbed particle-hole 
strengths are also shown in the lower panels for each nucleus. The arrows indicate the threshold 
energies for neutron emission. There are no significant strengths for K'^ = 3“ excitations in 
this low-energy region, so that they are not shown. Some particle-hole excitations near the Fermi 
surface are drawn by arrows with their values in the rightmost part for each nucleus. The solid, 
dashed, and dotted lines indicate the positive-parity levels, negative-parity levels, and the Fermi 
surface, respectively. The asymptotic quantum numbers [A^n^An] are indicated for pertinent 
levels. For the positive energy region, only the discretized continuum levels possessing resonance 
character are drawn. 
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Figure 7: Comparison of the SHF-RPA calculations with the Sill, SkM* and SLy4 interactions 
for low-frequency K'^ = 1~ and 2“ excitations on the SD states in 36,48,50g_ SLy4, the SD 
local minimum does not appear in Only excitations having greater than 10 W.u. 

are displayed here (see Fig. 6 for other excitations having smaller strengths). The numbers on 
the right-hand sides of individual levels indicate the K quantum numbers. The numbers beside 
the arrows indicate the intrinsic transition matrix elements squared of the isoscalar octupole 
operator, B{Q^^3), in units of fm®. Those for the electric octupole operator, B{E3), are also 
given in brackets in units of e^fm®. 
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